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Abstract 

We present analytical results for the 0(ot 2 s ) contributions to the functions t)a and rjv which 
parameterize QCD corrections to semileptonic b — > c transitions at zero recoil. Previously 
obtained approximate results are confirmed. The methods of computing the relevant two- 
loop diagrams with two mass scales are discussed in some detail. 
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1 Introduction 



The studies of the semileptonic decays of the 6-quark provide the best opportunity to de- 
termine | Kb |, a parameter of the CKM matrix and a fundamental input parameter in the 
Standard Model. Currently, two different methods of extracting the \V c b\ from the experi- 
mental data are used. One of them is based on the exclusive decays B — > D*lv\. In this 
method the recoil spectrum of the B meson decay is measured: 

— ^— = f(m B ,m D *,w)\V cb \ 2 F 2 (w) 1 + - In— ^ . (1) 

CLW V 7T TUB J 

In this equation w is the product of the four velocities of the B and D* mesons and / is a 
known function of the masses of the observed particles. 

In the limit when the b- and c- quark masses are considered to be large in comparison with 
Aqcd, the Heavy Quark Effective Theory |1], H, @J can be used to determine the nonper- 
turbative corrections responsible for deviations of .F(l) from unity in a model independent 
way. For this purpose, small velocity sum rules provide a reliable theoretical framework || . 

From the experimental point of view, the zero recoil point is not directly accessible due to 
the vanishing phase space. Therefore, one measures \V c b\ 2 ^(w) for w^l, extrapolates it 
to w — 1 and uses theoretical predictions for .F(l) to extract the value of \V c b\. 

Apart from the abovementioned nonperturbative effects, J-"(l) receives also perturbative 
QCD corrections. The Lorentz structure of the b — > c decay vertex is = 7^(1 — 75). We 
describe the modifications of the axial and vector parts at zero recoil by two functions: 

7 M -> Win, 7^75 -> ?M7/,75- (2) 
Both functions i]a and t]y can be expanded in power series in the strong coupling constant: 
VA,v = 1 + -C FV % + (-) 2 C FV % + O(o2). (3) 

The 0(a s ) effects have been known for a long time || [7], fj|. The 0(af) corrections to zero 
recoil form factors have been subject of controversy over several years. The need to perform 
complete two-loop calculation has been emphasized by many authors. The problem was 
solved in ref. [§], where 0(a 2 s ) were calculated in the form of a series in 5 = 1 -m c /m(,. The 
value of 0(a 2 s ) turned out to be relatively small and therefore the important source of the 
theoretical uncertainty in the |V^,| determination from exclusive decays has been removed. 

The purpose of the present paper is to present analytical results for the 0(a 2 ) corrections to 
two functions r] A and rj V and to discus the methods which have been used for their evaluation. 
The Feynman diagrams which have to be computed are shown if Fig. 1. 

There are several motivations for taking up this issue again. First we want to demonstrate 
that with the current techniques for doing multiloop calculations the calculation of 0(a 2 ) 
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(d) (e) (f) 

Figure 1: Two-loop QCD corrections to the b — > c transitions at zero recoil. Symbols <8> 
mark places where a virtual W boson can possibly couple to the quark line. 



correction to t]a,v is feasible and even relatively simple, if organized properly. Second, having 
a complete formula one can study it in the whole range of m c /mb, while the previously 
available expansion is not valid for m c <C m^. Finally, it is worth reminding that there are 
almost no results on the decays of fermions neither in QED nor in QCD at the two loop level. 
As a long standing example, the two loop QED corrections to the muon decay are not known. 
Being a challenging theoretical problem, the value of 0(a 2 ) corrections to muon decay does 
not influence precision electroweak physics at present. However, this situation may change 
if the electroweak parameters, such as sin 2 9w, are measured with higher accuracy, e.g. at 
the future linear collider. In view of the anticipated progress in precision measurements, 
accumulating as much information as possible about higher order effects in the fermion 
decays is warranted. 

The feasibility of the analytical calculation at zero recoil is mainly connected with the par- 
ticular kinematical configuration. We denote the momentum of the 6-quark in the initial 
state by p with the on-shell condition p 2 = m 2 . At zero recoil, the c-quark in the final 
state has the momentum up, where u> = m c /mb. The virtual W boson, which induces the 
semileptonic transition, carries off the momentum {l — u)p. Hence, the zero recoil kinematics 
dictates that the momenta of the particles in the final state are parallel to the momentum 
of the particle in the initial state. 

To put this in more simple terms, let us consider the zero recoil b — > c transition in the rest 
frame of the b quark. In this case the c-quark in the final state remains at rest and the 
energy release rrib — m c is carried off by the current which induces this transition. From this 
picture it is evident that no real gluon radiation can appear in the zero recoil transition. 

From the technical point of view, the fact that all momenta are parallel to each other 
simplifies significantly the calculations of the Feynman integrals. The Feynman integrals 
which appear are similar to those studied in the context of the quark wave function and mass 
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renormalization constants calculation, except that here we have to deal with two (rather than 
one) mass scales. Below we shall demonstrate that if this similarity is properly exploited, 
the analytical calculation of the form factors at zero recoil turns out to be feasible and even 
relatively simple. 

The paper is organized as follows. In section ^ we describe in some detail methods of 
computing two-loop integrals with two mass scales needed in the zero recoil calculation. We 
pay special attention to the subtraction of infrared divergences so that the difficult parts 
of the calculation can be done in 4 dimensions. In section |3] the special case of diagrams 
with a massive fermion loop inserted in a gluon propagator is considered. Our final result is 
presented and compared with the previous approximations in section £|. 



2 Basic scalar integrals 

To extract a contribution of a given Feynman diagram to rjAy we multiply it by u^^u^ 
or - u c 7^75Mfe and average over fermion polarizations. In this way one gets scalar products 
of four momenta in the numerator which one expresses through the available denominator 
structures. This leads to a reduction of the number of the propagators for a given Feynman 
diagram. After this step the remaining Feynman integrals have to be calculated. 

Before demonstrating various types of integrals we introduce some notations. We denote 
the loop momenta by k and I and work in the Euclidean space. There are nine possible 
propagators which appear in the Feynman integrals is this calculation: 

S 1 = k 2 , S 2 = l\ S 3 = (k-l) 2 , (4) 

S A = k 2 + 2pk, S 5 = l 2 + 2pl, S 6 = {k + l) 2 + 2p{k + l), 
S 7 = k 2 + 2upk, S 8 = l 2 + 2tupl, S 9 = (k + I) 2 + 2ujp(k + I). 

An arbitrary Feynman integral constructed from these objects will be denoted as 

D{a 1 , a 2 , a 3 , a 4 , a 5 , a 6 , a 7 , a 8 , a 9 \p, u) = / (5) 

J D 1 D 2 «J 3 4 D 5 Oq D 7 Dg Dg 



with 



[dk] = D = 4 - 26. (6) 



7T 



If the last two arguments of D(....\Q,x) are shown explicitly they denote objects which 
should replace p and u in Otherwise it is assumed that Q = p and x = u>. We use m h as 
a unit of mass. The proper dimension of the integrals can be easily restored, if needed. 
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For the presentation of our results for the relevant Feynman integrals it is convenient to 
introduce two auxiliary functions: 

7T 2 

R-y = — - Li 2 (l -uj) = Li 2 (cu) + ln(l - a;) ln(cu), 

R 2 = Li 2 (-w)+ln(l + w)lii(w). (7) 

If the argument of i?i j2 is given explicitly, it replaces u> in the above definitions. 



2.1 Scalar integrals with six propagators 

The highest number of propagators in Feynman integrals in the zero recoil calculation is 
six. Both planar (like in Fig. la,b) and non-planar (Fig. lc) integrals with six propagators 
appear. 



Planar diagrams with six propagators 

Let us consider the integral D{\, 0, 1, 1, 1, 0, 1, 1, 0): 



D(l, 0,1, 1,1, 0,1, 1,0) = J 



[dk][dl] 



k 2 (I - k) 2 (k 2 + 2pk)(l 2 + 2pl)(k 2 + 2upk)(l 2 + 2(vpl)' 



(8) 



The idea which permits the calculation of this integral is the following. Consider the product 
of 64 and S 7 or of S5 and S 8 in the integrand. For both it is possible to perform partial 
fractioning, like 



1 



1 1 



1 



(k 2 + 2pk) (k 2 + 2upk) 1 - uo k 2 \ k 2 + 2pk k 2 + 2uopk J ' 
For D(l, 0, 1, 1, 1, 0, 1, 1, 0) we therefore write 



(9) 



S4S7S5S8 (1 — lo) 2 S1S2 \Sa Sj J \Ss S; 



(10) 



and, accordingly 



D(l, 0,1, 1,1, 0,1, 1,0) = 



1 



D(2, 1, 1, 1, 1, 0, 0, 0, 0) + u 2 D(2, 1,1,0, 0, 0, 1, 1, 0) 



(1-u) 2 

-uD(2, 1, 1, 1, 0, 0, 0, 1, 0) - cuD(2, 1, 1, 0, 1, 0, 1, 0, 0) 



(11) 



We notice that the first two D structures in the square brackets are two-loop two-point 
on-shell integrals with a single mass scale. The only difference between them is in the 
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incoming momenta and in the masses of the particles inside the loop. An algorithm to 
compute an arbitrary integral of this type is known |], [IJj, [IT] (the underlying idea is the 
integration by parts method ||12|| ). A large number of such integrals had to be computed in 
the approximate approach to the present problem || . For that purpose a new implementation 
of the recurrence algorithm was necessary. The ability to compute such integrals determines 
our strategy also in the present case. 

Apart from the known single scale integrals, eq. (|TTJ) contains two unknown integrals. Let 
us consider one of them as an example: 



D(2, 1,1, 1,0, 0,0, 1,0) = J 



[dk)[dl] 



(k 2 ) 2 l 2 (l - k) 2 (k 2 + 2pk){l 2 + 2upl) 



The price we have to pay for partial fractioning is that the integrals we obtain have stronger 
infrared divergences than the original ones. The integration in D-dimensions is much more 
complicated than if the number of space-time dimensions could be equated to 4 from the 
start. In circumventing the problem the following trick proved to be useful. We introduce a 
Feynman parameter to combine two propagators in the integrand of-D(2,l,l,l,0,0,0,l,0): 

' dt (1 - t) 



(k 2 ) 2 (k 2 + 2pk) J (k 2 + 2pkt) 3 
Introducing the variable t = ufi one gets 

dt (1- t) 1 f J du(l - uu) 



(k 2 + 2pktf J (k 2 + 2pkujnf 

Splitting the integration over /i into two parts: from to 1 and from 1 to 1/w, we get the 
following representation for D(2, 1, 1, 1, 0, 0, 0, 1, 0): 

D(2, 1, 1, 1, 0, 0, 0, 1, 0) = uD(2, 1, 1, 0, 0, 0, 1, 1, 0) + cu(l - u)D(l, 1, 1, 0, 0, 0, 2, 1, 0) 

+w(l - cu) 2 D(0, 1, 1, 0, 0, 0, 3, 1, 0) + A, (12) 



where A is: 

rWi n(o i 1 n i n 3 n nl™,» ,A — 



A = 2u J dfi(l - ufi) 



D(0, 1, 1, 0, 1, 0, 3, 0, 0\pu, n) - ^D(0, 1, 1, 0, 1, 0, 3, 0, 0\pu, 1) 



The advantage of this representation is that all the D functions, listed in eq. (12) are now 
of the single scale type and their values are available. On the other hand, the part of the 
result denoted by A is finite and its calculation can be performed in four dimensions. 

Let us stress this point once again. One of the obstacles in calculating the zero recoil form 
factors is the fact that individual Feynman integrals have infrared divergences, which should 
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disappear in the sum. It is the appearance of these divergences, which makes the calculation 
tough. We have demonstrated above that in the zero recoil calculation one can extract the 
infrared divergences in the form of the on-shell integrals which are known. 

For the sake of illustration, we write down the following finite result: 

D(0, 1, 1, 0, 1, 0, 3, 0, 0\pu, fi) - ^D(0, 1, 1, 0, 1, 0, 3, 0, 0\pu, 1) = ^i^Z^f- ( 13 ) 

The simplicity of this result leaves no doubts that the remaining integration over fi can easily 
be performed. We thus arrive at the result for the planar diagram with six propagators: 

D(l, 1,0,1, 1,0, 1,1,0) 
1 



to 2 



12e 2 + 2e-l l + u T ., N (2 + u) \n 2 (u) (1 - 2e) ln(w) 13vr 2 

Ll 2 (l - UJ) =7- ; h 



k 2 1-u K ' 2(1-0;) 26(1-0;) 48 



(14) 



This is the only planar diagram with six different propagators which is necessary for the 
present calculation. 

Non-planar diagrams with six propagators 

Another integral with 6 propagators we need is D(l, 1, 0, 1, 1, 1, 0, 1, 0). Performing partial 
fractioning we get 

D(l, 1, 0, 1, 1, 1, 0, 1, 0) = — !— 2, 0, 1, 1, 1, 0, 0, 0) - j^—Dil, 2, 0, 1, 0, 1, 0, 1, 0). 

1 — CO 1 — UJ 

The first of these integrals is of the single scale type. To calculate the second one, we use 
the same trick as described above and get the following representation: 

D(l, 2, 0, 1, 0, 1, 0, 1, 0) = -D(l, 2, 0, 1, 1, 1, 0, 0, 0) 

UJ 

-i^D(l, 1, 0, 1, 2, 1, 0, 0, 0) - - / ^ ( 1 - -) D(0, 1, 0, 3, 0, 0, 0, 1, l]*" 1 ). (15) 

LU Z Ul J t \ UJ J 

LO 

The first two integrals in this expression are of the single scale type. The expression for the 
Feynman diagram under the t integration will be presented later (eq. |2^). Direct integration 
gives: 

1 1 (i _ 1) m x, o, 3, o, o, o, i, i|r>) = -^^ft + 

1 - <-. _ A laH _ ZOzii _ AizM. (16) 



4a; 8uj 8a; 48a; 



From this an expression for D(l, 1, 0, 1, 1, 1, 0, 1, 0) can be constructed. 

The diagram .0(1,1,0,1,0,1,0, 1,1) represents a more challenging problem. We note that 
it has an overall infrared divergence and cannot be calculated in D = 4. To compute it we 
introduce the operator 



d 



For dimensional reasons the integral we are interested in is proportional to (p 2 ) 2 2e . By 
applying the operator L to it one finds 

LD(1, 1, 0, 1, 0, 1, 0, 1, 1) = -(4 + 4e) D(l, 1, 0, 1, 0, 1, 0, 1, 1). (17) 

By applying the same operator to the integrand we get 

- (1 + At) D(l, 1, 0, 1, 0, 1, 0, 1, 1) = D(0, 1, 0, 2, 0, 1, 0, 1, 1) + D(l, 0, 0, 1, 0, 1, 0, 2, 1) 



1-uj 



D(l, 1,0, 1,0, 2, 0,1,0) 



UJ 



1-uj 



D(l, 1,0, 1,0, 0,0, 1,2). 



(18) 



It is useful to note that Feynman integrals in the above equation satisfy the following sym- 
metry relations 

D(l, 0, 0, 1, 0, 1, 0, 2, 1) = Lu~ 4 ~ 4e D(0, 1, 0, 2, 0, 1, 0, 1, 

D(l, 1, 0, 1, 0, 0, 0, 1, 2) = cj~ 4 ~ 4e D(l, 1, 0, 1, 0, 2, 0, 1, Ol^r 1 ). 
Therefore, it is enough to know two of them. 

We present the calculation of D(0, 1, 0, 2, 0, 1, 0, 1, 1) which turns out to be the most compli- 
cated integral. Applying the operator L to it once more, we get the following identity: 



-2£(0, 1,0, 2, 0,1, 0,1,1) = 



1-uj 



D(0, 1,0, 3, 0,1, 0,1,0) 



-— ^— O>(0, 1, 0, 3, 0, 0, 0, 1, 1) + -D(0, 1, 0, 3, 0, 1, 0, 0, 1) - -D(0, 0, 0, 3, 0, 1, 0, 1, 1) 

1 — UJ UJ UJ 

+ -^—D(0, 1, 0, 2, 0, 2, 0, 1, 0) + D(0, 0, 0, 2, 0, 1, 0, 2, 1) - -^O>(0, 1, 0, 2, 0, 0, 0, 1, 2). 
1 — uj 1 — UJ 

(19) 



All Feynman integrals in this equation are finite and relatively simple to calculate directly. 
The results for individual Feynman integrals read: 



D(0, 1,0, 3, 0,0, 0,1,1) = — 



£(0,1,0,2,0,0,0,1,2) 



UJ' 



4^ 2 _ Rl ) + 2(T^) 



+ 



71 

+ 



(1-LU) 2 2(1 -d) 8 

u;ln(cj) 1 7r 2 



(20) 



+ 



{1-uj) 2 2(1 -u) 8 
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D(0, 1,0, 3, 0,1, 0,0,1) = 
D(0, 0,0, 3, 0,1, 0,1,1) = 
D(0, 0,0, 2, 0,1, 0,2,1) = 



1 — uj + \n(uj) 



2(1 -uj) 2 

1 + UJ 



uj ln(u;) 



7T UJ 



uj(1 — uj) 

1+OJ 

~uj 2 {1-uj) 



8^ Rl R ^> 4(1 + uj) 2 161 + a; 



uj\yl(uj) a; 
?(l + cu) 2 + 



8 1 + uj 



2D(0, 1, 0, 3, 0, 1, 0, 1, 0) + D(0, 1, 0, 2, 0, 2, 0, 1, 0) = -j-^ - i? 2 ). 

The last missing contribution for the scalar non-planar diagram with six propagators is 
D(l, 1, 0, 1, 0, 2, 0, 1, 0). This integral is infrared divergent and therefore its direct evaluation 
is not convenient. However, the trick described in the context of the planar integral with six 
propagators is useful also here. One finds the following representation: 



D{1, 1,0, 1,0, 2, 0,1,0) = - 

UJ 



r dt 

D{\, 1, 0, 1, 1, 2, 0, 0, 0) - J -D(0, 1, 0, 2, 0, 0, 0, 1, 2\r l ) 



The result for the Feynman integral under the integration over t is presented in the the list 
of the integrals given above. The final integration over t is elementary. As a result we find 



D(l, 1,0, 1,0, 1,0, 1,1) 



+ 



4oj 2 e 4uj 2 (1-uj) 



2(1 + u) - R 2 ) + 4 ln(uj) + 12(1 - uj) - n 2 



[21) 



2.2 Scalar integrals with five propagators 

The scalar integrals with five propagators represent probably the most difficult part of this 
work. The reason is that these integrals contain trilogarithms, therefore evaluation of the 
corresponding Feynman integrals is harder and the resulting expressions are rather compli- 
cated. 



Planar diagrams with five propagators 

The most difficult planar diagram is D(l, 1, 1, 1, 0, 0, 0, 1, 0). To calculate it we use the trick 
described in the previous section, which allows us to extract on-shell divergent piece from 
this diagram. The following representation is useful: 

D(l, 1,1, 1,0, 0,0, 1,0) = - 

UJ 



£>(1, 1,1, 1,1, 0,0, 0,0) -J dt D(l, 0,1,1,0,0,0,2,0|0 
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where 

2 , N ln 2 (t) 2tt 2 
D(l, 0,1,1, 0, 0, 0, 2, 0|t) = --i?i(t) - j±j + — . 

Performing the integration over t one gets: 

1 , x 6(Li 3 (w) - Ca) 
£1,1,1,1,0,0,0,1,0 = -D 1,1,1,1,1,0,0,0,0 ^ 

OJ 

1 



+ ■ 



2vr 2 

41n(w)Li 2 (cj) + ln 2 (u;) ln(l - u) + — \n 2 (oj) 

o 



(22) 



Non-planar diagrams with five propagators 



This is the most difficult case. Let us start our discussion with D(l, 1, 0, 1, 0, 0, 0, 1, 1). We 
find that this diagram satisfies the following relation: 

L>(1,1,0,1,0,0,0,1,1) = ^#(-), (23) 

OJ 2 \ujj 

where the function H is given by 

1 ^ 

H(u) = - f D(l,0,0,l,0,l,0,2,0|*)d*. 



For the expression under the integral we get the following result: 

D(l, 0, 0, 1, 0, 1, 0, 2, 0|f) = - 2 ^±^R 2 (t) + ^ - (24) 

This nice integral representation essentially solves the problem. The subsequent integration 
over t can be performed; however, the final result is complicated and we relegate it to the 
appendix. For all practical purposes, like evaluation of the final result as a function of oj, 
the integral representation provides a more useful starting point than the exact formula. 

There are two other integrals which are necessary. We will discuss the calculation of one of 
them D(0, 1, 0, 1, 0, 1, 0, 1, 1). For the time being let us denote it as F(oj). 

It turns out to be convenient to derive a differential equation for this function. Differentiating 
it with respect to oj one gets: 

; d ' ' " LJ F(oj) = J(oj), (25) 



\ doj oj(1 — oj) 

J{oj) = --D(0, 0, 0, 1, 0, 1, 0, 2, 1) + — !—£>((), 1, 0, 1, 0, 0, 0, 1, 2). 
oj 1 — OJ 
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Writing F(uj) = p(uS)G{uj) with 

p{u) 



uj(1 — uj) 
one gets a new equation for G(uj): 

&G(uj) 



duj 



u(l-u)J(u). (26) 



In order to determine the right hand side of this equation we have to find J(uj). For this we 
need two Feynman integrals which can be computed directly 

D(0, 1,0, 1,0, 0,0, 1,2) = _^l + 4L±4 fl2 _?^!M + ^-) 



uj uj(1 — uj) 1 — UJ 3uj(1 — uj) 

D (o,o,o,i,o,i,o,2,i) = ?ft_^ + (i + Mi.'M _^) 

uj 1 — uj 1 — uj 2 3(1 — uj 2 ) 

and 



2R X AR 2 (l + 3u)\n 2 (u) 2it 2 (2-uj) 
1 — 1 — uj z 3(1 — ur' 



With this input the differential equation (|26|) leads to an integral representation for F(uj). 



i 

D(0, 1,0, 1,0, 1,0, 1,1)= 1 / dtJ x (t). 

uj(1 — UJ) J 



uj(1 — uj) 

The analytical result can be found in the appendix. 

2.3 Feynman integrals with four propagators 

A lot of integrals with four propagators appear in the course of this calculation. A large 
fraction of them can be calculated directly. In some cases however some special methods 
appear to be useful. 

Below we list some useful results. For the diagram D(0, 0, 0, 1, 0, 1, 0, 1, 1) we find the fol- 
lowing: 

w o, o, i, o, i, o, i, i) = (1 -»'^ r <f>f 7 1 -'* - i-L- U M - aFl a 

e (1 — 6e) (1 — uj) 1 — uj \ \uj 

where 

„, . 2(1 -uj) 2 R 1 2(1 -uj)R 2 uJ 2 \n 2 (uj) ti 2 (1 + 2uj 2 ) 
Pl(uJ) = h 



uj uj 1 + uj 3(1 + uj) 
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For the diagram D(l, 0, 0, 1, 0, 1, 0, 1, 0) we have: 

r(2e)5(l-e,l-e) 



D(l, 0,0, 1,0, 1,0, 1,0) 



e 1 -3e 



2(1 + U ) Ri _ 2(1 + uf R ^ _ u; 2 ln 2 M _ tt^I+oO _ (2f 



LO LO " 1 — UJ 



Another type of diagrams we have to deal with are those with an internal massive bubble, 
for example D(0, 1, 0, 1, 0, 1, 0, 1, 0). Similar integrals were studied in |J in connection with 
the investigation of the quark wave function and mass renormalization constants. Methods 
developed there permit a calculation of the integrals of this type necessary in our case and 
will be useful also in case of massive quark loop insertion in the gluon propagator (presented 
in the next section). For example for the diagram D(0, 1, 0, 1, 0, 1, 0, 1, 0) we find 

D(0, 1,0, 1,0, 1,0, 1,0) = £>(1, 1, 1,0, 0,0, 0,1,0) + 2(1 Rl 

LO 1 6 



Other diagrams which appear in the course of this calculation, like the diagrams with three 
propagators, can be reduced to the diagrams with four propagators similar to the ones 
discussed above but with the different structures in the numerator. All of these diagrams 
can be studied in a similar manner and we do not present them here in any detail. 

The results presented above are sufficient to get the analytical result for the major part of 
the diagrams involved in this calculation. However, they are not sufficient for the diagrams 
with the internal massive fermion loop. We shall discuss their evaluation in the next section. 



3 Diagrams with a massive fermion loop 

The diagrams in Fig. [TJ with a massive b or c quark loop insertion in the gluon propaga- 
tor represent a somewhat special case. It is clear how to extract their divergences. The 
ultraviolet divergences are canceled by wave function renormalization. On the other hand, 
these diagrams give infrared finite contributions if we use on-shell renormalization for the 
coupling constant. Therefore, in order to deal with the finite contributions from the very 
beginning one should study the correction to the vertex together with the wave function 
renormalization and renormalize the coupling constant in a QED-like way. 

The contribution of the massive fermions to the gluon polarization function U(k 2 ) can be 
written in a form of a dispersion integral (for the purpose of illustration we consider a particle 
in the internal loop to be 6-quark) subtracted at k 2 = 0. After that, the result appears as 
a convolution of a one-loop correction evaluated with the gluon of mass A with the fermion 
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spectral density. The following transformation rule is valid: 



d\ 2 



1 + 



2ml 



4mr 



A 2 (P - A 2 ) \ A 



4m 2 
\2 • 



(30) 



To proceed further, we first calculate the one-loop form factors t]a,v with a massive gluon. 
Afterwards the result must be integrated over the mass of the gluon with the fermion spectral 
density: 



VA,v{b) 



d\ 2 



i + 



2m 2 b 



1- 



4m 



A 2 



6 (1,A) 

- Va,v 



(31) 



lint 



Here riA,v{b) denotes the contributions of diagrams with a b quark loop inserted in the gluon 
propagator and % v denotes the one-loop form factors calculated with a massive gluon 
propagator l/(k 2 — A 2 ). 

It is advantageous to write the one-loop result leaving the integration over the last of the 
Feynman parameters untouched. Namely, both one-loop form factors calculated with the 
finite mass of the gluon ry^v can be written as: 



Vav 



d(3\ 



where 



p^m 
p 2 a (p^)-- 



X 2 (3 + (l-/3) 2 \ 2 f3 + uj 2 (l- (3) 2 1 

lu(2 - 29/3 + hl(3uj + 34/3 2 - 30(3 2 lu - 7(3 3 + 9(3 3 lu + 18a;) 

6(1 -w) 

oo 3 (30(3 2 - 9(3 3 + 7(3 3 u - 34(3 2 uj -2u- 51/3 + 29pu - 18) 

6(1 -uj) 

uj{-2 -7(3 + 17(3uj + IA(3 2 - 1®(3 2 uj - h(3 3 + 3(3 3 u + 6u) 

2(1 -w) : 
cj 3 (-6 + 2u - 17(3 + 10(3 2 - 3(3 3 + 7(3u - U(3 2 u + 5(3 3 uj) 



(32) 



2(1- 

Substituting this result to eq. (|3lD and performing the integration over A one gets 



VA,v(b) ~ d(3 



Ay i 



n 







uj 2 {1-(3) 2 



n 







where II (s) is the vacuum polarization function: 



U(s) 



1 ) v In 

s 



r : 1 x 4 5 



v — 1/ s 3 



(33) 
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with 

Subsequent calculation of the one-parametric integrals is "tedious but straightforward." 

So far we have discussed the calculation of the massive fermion bubbles in what can be called 
on-shell or QED like normalization. The final results which we would like to get are the ones 
with the MS'-coupling constant a s renormalized at the symmetric point /x = ^Jmtfn~ c . The 
transformation to this form amounts to performing a scale transformation in the coupling 
constant. Finally, for the tVquark bubble contribution to the zero recoil form factors one 
gets: 

(3c; 4 - 7uj 3 - 7uj 2 - 25c; + 36)R 1 (5uj - 18c; 2 - 34c; 4 - 21c; 5 + 36)R 2 
VA[b) = T^— d + 



48c; 4 48o; 4 (l - uj) 

1 + u) \n 2 (uj) (15 - 4c; + 18c; 2 - 37c; 3 ) ln(o;) 72 - 33c; + 104c; 2 - 143c; 3 tt 2 (5 - 3c;) 



4(1 -w) 24(1 -c;)c; 2 72o; 2 (l - uj) 36(1 - u 

and 

(c; 4 + 3c; 3 + 3c; 2 - I9u + 12) R x (-7c; 5 + 12 - 22c; 2 - 9c; - 6cu 4 )R 2 



Vv(b) 



16c; 4 16c; 4 (l - uj) 



to 



l+c;)ln 2 (c;) (15 - 36c; + 14c; 2 - 25c; 3 ) ln(c;) 24 - 51c; + 62c; 2 - 35c; 3 tt 2 (3 
4(l-a;) 24c; 2 (1 - uj) + 24c; 2 (l - uj) 12(1 - uj) 



4 Results 



(2) 

It is convenient to divide up the functions ry A V into parts proportional to various SU(3) 
factors (an overall factor Cp has been factored out): 



(2) 

Va!v 



CfVa,v + (Ca 



1C F )r,i F v + T R N L -nh v + T R r,% 



\ ■ 



(34) 



For a general SU(N) group C A = N; C F = (N 2 - 1)/2N; T R = 1/2. N L denotes the 
number of the light quark flavors whose masses can be neglected. The last term contains 
contributions of the massive quark loops, with b and c quarks. We neglect the top quark; its 
impact is suppressed by a factor ~ m^/m 2 . We use the pole masses of b and c quarks and 
express our results in terms of ajjg(y/m b m c ). 

Among the eight coefficient functions in eq. fl34j) i]\ v are already known exactly |13|| . They 
corresponds to our diagram (f) in fig. 1 with a massless fermion in the loop; they read 



1 

24 
5 

24 



1+UJ 



1-uj 

1+UJ 



\n(uj) + 2 



44' 

ln(c;) H 

1-uj K 1 15. 



(35) 
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The contribution of the massive quark loops is obtained by combining the result for the b 
quark loop, obtained in the previous section, with an analogous formula for the c quark. 
The latter is obtained from the former by a simple symmetry transformation. We find 

H „ N3/1 , 12 + 5a> + 13c; 2 + 5c; 3 + 12c; 4 n 
Vv = (1 - ^) 3 (1 + ^ Ri 

, 12 - 21u - c; 2 + c; 3 - 14a; 4 + c; 5 - o> 6 - 21c; 7 + 12c; 8 n 

16(1-^ * 2 

,, x 1 + c; + 5c; 4 - 3c; 5 , , w x 15 - 51c; + 40c; 2 - 51c; 3 + 15c; 4 



+ ln » S ln(a;)(l+a;) , 

4(1 — c;) 24(1 — u;jur 

37 - 30c; + 66c; 3 - 53c; 4 + 12c; 5 , 24 - 27c; + 70c; 2 - 27c; 3 + 24c; 4 

-7r + oFT? > ( 36 ) 



96(1 -a;) 24c; 2 



„ , , _ 2l 36 - 25c; + 29c; 2 - 32c; 3 + 29c; 4 - 25c; 5 + 36c; 6 
' 48c^ 

36 - 31c; + 13c; 2 - 13c; 3 - 42c; 4 - 13c; 5 + 13c; 6 - 31c; 7 + 36c; 

48(1 -c;)c; 4 



Va = (1 - ^) ! ! Ri 

1 + a;) — — ' — — /7 ' ~~ ' — R2 



, 2/ N 3 + 3c; + 7c; 4 - 9c; 5 , . s 15 - 19c; - 19c; 3 + 15a; 4 

+ ln 2 (c;) ln(c;)(l + u 



12(1 -a;) " " 24(1 -c;)c; 2 

4 _L q«,.,5 



72 + 39c; + 286c; 2 + 39c; 3 + 72c; 4 79 - 10c; + 54c; 3 - 127c; 4 + 36c; 5 



-It 



(37) 



72c; 2 288(1 - u) 

The functions R\ 2 are defined in eq. (|7|). 

Diagrams without fermion loops give rise to two color structures and we chose the combina- 
tions Cp and CaCf — 2C\ to describe them. For the functions T)y A we find 

F a> 2 (2 + a>), 2 / \ /89 tt 2 \ 1 + . , (1 + c;) 3 T . „ N 
l-10u, + ^ v + 53 + 3-10u,H-3 W - 



24c; 48 (1 - ^ 

p c; 2 (6-c;) , 2 , , /53 tt 2 \1 + c; 1/n (1 + u)(l - 6c; + c; 2 ) T . ,„ 



1 - 38c; + c; 2 2 _ 95 9 - 14c; + 9a; 2 
72^ ^ ~72 + 3(1 - c;^ 2 



7r 2 -- + — -^—f(u), (39) 



7T 2 \ , , s 3 3 /, , , 1\. , 



/(a;) = ( Li 2 (l - w) - — 1 In(w) + - (Li 3 (c;) - Cs) + ^ ^(1 - u>) + - J hi 2 (c;). 

The most complicated and difficult to compute are the functions 77^4. We parameterize them 
using, in addition to the notations introduced before, in terms of four functions /i...4(u;) for 
which closed formulas are given in the appendix. In the following formulas we drop the 
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argument u of ff 

af q \ , n . A + 2uj-uj 2 1 - a; + 5a; 2 - a/ + 

=— {1 - 3uj) - h{1 + Uj) ^^ h 16(l-a;) +/4 80^) 

(1 - w 2 )(3i2! + ffc) _ ^ ^ + 3 _ + 3 _ ^ _ 1 ^ 



16a; v ' 2 v v y ^ ' 4 

17, , x l + w , . x (l + a;) 2 7r 2 5 . . , 17 

-- In a; In a; V ; + - 1 - a; tt 2 - -, 40 

96 1 — u 24 1 — a; a; 96 48 



r^ F = -R, \n(u) + ^-U 3 (u) + ^ /v ~ ~ ' ~~ ' Ri 



3 T . , , 1 - a; 2 8-a; + 8a; 2 ^ 

-Li 3 (a;) + =- 

2 3V ; 48a; 2 



(l+a;)(8-3a; + 22a; 2 -3a; 3 + 8a; 4 l 3, , . ftw, 

+ " 48(1 -„V + l' n(1 " W) k (W) " + M 

3 + u + 19a; 2 - 3a/ 3- lla> + 31a; 2 -3a; 3 a;(7 + 3a>) 

~ h 48^ h 48(1 -u) + /4 24(l-cu) 

, ,, , 9 - 6a; + 5a; 2 - 4a> 3 2l , N 3-2a; + 3a; 2 61 , , . 1 + u 

— In (a; la; : tt; vr Info;) ; r mlu) 

V ; 24(1 -a;) 2 V ; 72(1 - u)u 96 K ' 1 - u 

63 - 46a; + 23a> 2 - 8a; 3 „ 151 7r 2 ln2 5 A 

H ; ; vr 2 (3- (41) 

288(1 -a;) 72 6 4 s V ; 

The behavior of the above coefficient functions is plotted in Fig. |2|. It can be seen that the 
corrections to the vector coupling vanish when u — ► 1, in accordance with vector current 
conservation. For ui = 0.3, which roughly corresponds to the physical value of the mass 
ratios of c and b quarks, we find the following results 

r$ = 0.107667, 7$ = 0.393822, t]$ f = -0.167695, 

rfi = -0.154818, rf A = -0.586889, r^ F = -0.909398, (42) 

in excellent agreement with the approximate values found in || (see eq. (11) in that refer- 
ence). 

From the above exact formulas we can also find the expansions of all the coefficient functions 
around ui = which corresponds to the vanishing c quark mass. We give below the terms 
which do not vanish as u — > 0: 



ln 2 (o;) 131n(o;) 1147 37vr 2 



ji 



4 24 ' 288 96 ' 

ln 2 M 7 Info;) 971 79 tt 2 



4 8 288 288 

f 91nV) /113 tt 2 \ 41 7tt 2 9Cs 

* -^2— + (^-T2j lnH + 48 + l2--— ' 

IA 32 V 96 12 / v 7 72 36 2 
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(a) (b) 

Figure 2: Coefficient functions rf (solid), r\ H (dashed), r\ AF (dotted), for the vector (a) and 
axial case (b), plotted against uj = 0..0.8; none of these functions varies noticably between 
oj = 0.8 and lu — 1. 



5 


n 2 \ 


ln(w) 


23 


TT 2 


96 " 


^ 24 J 


~ 48 


96 4 


49 




InH 


83 


5tt 2 
+ 32 


96 " 


^ 24 j 


~ 36 



vr 2 ln(2) 9C 3 



2 4 



^ - (-^ + ^)lnH-^ + ^r + ^^-2C3. (43) 



5 Conclusions 

In this paper we have presented an analytical calculation of the 0(a 2 s ) corrections to the 
functions 7]a and r]y which parameterize semileptonic b — > c transition at zero recoil. Nu- 
merically, our results are in excellent agreement with the results of the ref. H obtained using 
an approximate approach to the present problem. 

This calculation was facilitated by an efficient method of computing two-loop vertex functions 
with two mass scales in the zero recoil limit which we demonstrated in some detail. Infrared 
singularities of individual diagrams were shown to be the same as in certain combinations 
of single scale integrals, for which an exact algorithm is known. We presented various useful 
tricks which permit a calculation of the remaining two-scale integrals. Since the previous 
approximate formulas H were valid only for similar masses of initial and final quarks, the 
exact formulas we obtained significantly extend our knowledge of two-loop corrections to 
fermion decays. 
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A Analytical results for the functions fi 



In eq. (f40l - f41~|) we presented the results for the coefficients rj AF in terms of four functions /1...4. 
In this appendix we present closed formulas for f\ and f^. The remaining two functions f'2,3 
are obtained from these results by applying the following symmetry relations: 



UJ \0J 



UJ" \UJ 



For f x {u) = D(0, 1, 0, 1, 0, 1, 0, 1, 1) we find 
1 



-4Li? 



uj(1 — oj) 
1+oj 



4Li 3 (cu) + 8Li 3 (y^j) + 4Li 3(! - w) + 6Li 3 (-cj) - 2Li 3 (l - uj 2 ) 



\- I Li 3 ( tt—\ - ^Cs - 4Li 2 ( ) ln(w) - 6 ln(w)Li 2 (-w) - 2 ln(u;)Li 2 (u;) 



2 / Vi + 2 " VI + cj 

+2 ln(2) ln 2 (l + uj) - 2 ln 2 (2) ln(l + w) - 2 ln 2 (u;) ln(l - w) - \n 2 (uj) ln(l + w) 

7T 2 2 

+7T 2 ln(l + u) ln(2) + - ln 3 (2) - 2 ln 3 (l + u) 

3 3 



The result for f 4 (u) = D(l, 1, 0, 1, 0, 1, 0, 1, 0) is 



1 

UJ 



- 2Li 3 (l - cj 2 ) - 4Li 3 ( + 4Li 3 (l - + 8Li ' 1 



1 +UJ 



2U 3 (uj) 



+4Li 3 (^-) - ^-(3 - 4Li 2 (r^-) ln(w) + 2 ln(cu)Li 2 (cu) - 2 ln(w)Li 2 (-w) 

7T 2 2 

+2 ln(2) ln 2 (l + w) ln(2) + 2vr 2 ln(l + uj) - \ii 2 (uj) ln(l - oj) + - ln 3 (2) 

3 3 



-21n 3 (l + w) -21n 2 (2) ln(l + w) 
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